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With this work we elaborate on the physics of quantum noise in thermal equilibrium and in station- 
ary non-equilibrium. Starting out from the celebrated quantum fluctuation-dissipation theorem 
we discuss some important consequences that must hold for open, dissipative quantum systems 
in thermal equilibrium. The issue of quantum dissipation is exemplified with the fundamental 
problem of a damped harmonic quantum oscillator. The role of quantum fluctuations is discussed 
in the context of both, the nonlinear generalized quantum Langevin equation and the path inte- 
gral approach. We discuss the consequences of the time-reversal symmetry for an open dissipative 
quantum dynamics and, furthermore, point to a series of subtleties and possible pitfalls. The path 
integral methodology is applied to the decay of metastable states assisted by quantum Brownian 
noise. 



This work deals with the description of quan- 
tum Brownian motion in Hnear and nonUnear 
quantum systems that exhibit frictional influ- 
ences. The symmetries of thermal equilibrium 
impose severe constraints on the time evolu- 
tion properties of open quantum systems. These 
lead to a quantum generalization of the classical 
Einstein relation that connects friction with the 
strength of thermal quantum fluctuations. There 
exist a variety of theoretical roadways to model 
quantum dissipation. Here, we discuss the topic 
for the prominent case of a damped harmonic 
oscillator upon combining thermodynamics and 
linear response theory. A dissipative nonlinear 
quantum dynamics can be dealt with a general- 
ized quantum Langevin equation, a path integral 
formulation, or in terms of a generalized quan- 
tum master equation for the corresponding re- 
duced dynamics. We illustrate the situation for 
the problem of the dissipative decay out of a 
metastable state. Furthermore, we point out a 
series of subtleties, pitfalls and shortcomings that 
one must be aware of when confronted with the 
world of quantum noise driven phenomena. 



I. INTRODUCTION 

Albert Einstein explained the phenomenon of Brow- 
nian motion in his annus mirabilis of 1905 by use of 
statistical methods which he ingeniously combined with 
the laws of thermodynamicsi. In this pioneering work 
he as well provided a first link between the dissipative 
forces and the impeding thermal fluctuations, known as 
the Einstein relation which relates the strength of dif- 
fusion to the friction. This intimate connection between 
dissipation and related fluctuations was put on a firm ba- 
sis much later when Nyquisli^ and Johnson'^ considered 
the spectral density of voltage- and current-fluctuations. 

What role do quantum mechanics and the associated 
quantum fluctuations play in this context? After the 
birth of quantum mechanics in the early 1920's we can 



encounter in the very final paragraph of the 1928 pa- 
per by Nyquist for the first time the introduction of 
quantum mechanical noise via the substitution of the 
energy kT from the classical equipartition law^ by the 
thermally averaged quantum energy (but leaving out the 
zero point energy contribution) of the harmonic oscilla- 
tor. Nyquist's remark thus constitutes a precursor of the 
celebrated work by Callen and Welton^ who generalized 
the relations by Einstein, Nyquist and Johnson to include 
quantum effects; In their work they put forward a gener- 
ally valid connection between the response function and 
the associated quantum fluctuations in equilibrium, the 
quantum fluctuation-dissipation theorem. 

Without doubt, quantum fluctuations constitute a 
prominent noise source in many nano-scale and biolog- 
ical systems. For example, the tunnelling and the trans- 
fer of electrons, quasi-particles, and alike, is assisted by 
noise for which the quantum nature cannot be neglected. 
The features of this noise change drastically as a func- 
tion of temperature: At sufficiently high temperatures a 
crossover does occur to classical Johnson-Nyquist noise. 

With this work we shall present various methods and 
schemes of modelling quantum Brownian motion from 
first principles. In particular, the thermal noise must at 
all times obey the quantum version of the fluctuation- 
dissipation theorem (a la Callen- Welton) . This latter 
property is necessary in order to be consistent with 
the second law of thermodynamics and the principle of 
(quantum) detailed balance. We elaborate on several 
alternative but equivalent methods to describe quan- 
tum noise and quantum Brownian motion per se: These 
are the functional integral method for dissipative quan- 
tum systemSi^ and time-dependent driven quantum 
systemsS, the quantum Langevin (operator) approach^-, 
stochastic schemesiStii, or the concept of stochastic 
Schrodinger equations^^. In doing so, we call attention to 
distinct differences to the classical situation and, as well, 
identify a series of delicate pitfalls which must be ob- 
served when making even naive looking approximations. 
Such pitfalls involve, among others, the rotating wave ap- 
proximation, the use of quasi-classical Langevin forces, 
the quantum regression hypothesis and/or the Markov 
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approximatio: 
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II. THE QUANTUM FLUCTUATION-DISSIPATION 
THEOREM AND ITS IMPLICATIONS 

As already mentioned, in 1951 Callen and Welton 
proved a pivotal relation between equilibrium fluctua- 
tions and dissipative transport coefficients. Note also 
that this quantum fluctuation-dissipation relation holds 
true independent of particle statistics. The following 
cornerstone achievements can be found in this primary 
workS^: 

• The generalization of the classical Nyquist's for- 
mula to the quantum case. 

• The quantum mechanical proof that susceptibilities 
are related to the spectral densities of symmetrized 
correlation functions. 

For a single degree of freedom, linear response the- 
ory yields for the change of the expectation value of an 
operator- valued observable B due to the action of a (clas- 
sical) force F{t) that couples to the conjugate dynamical 
operator A 



<5B{t)> 



dsxBA{t - s)F{s) . 



(1) 



Here, 5B{t) = Bit) — <B>o denotes the difference with 
respect to the thermal equilibrium average <B>q in the 
absence of the force. The reaction of the system is con- 
tained in the response function XBA{t) with a so-called 
dissipative part 



XAB{-t)] . 



(2) 



The Fourier transform of XbaW '^^^1 be denoted by 
X^^(ijj). It is worth noting here that only when A — B 
does this part in fact coincide with the imaginary part of 
the complex- valued susceptibility x'^^(aj). 

The fluctuations are described by the equilibrium cor- 
relation function 



CsAit) = <SB{t)SA{0)>,3 



(3) 



at inverse temperature P — 1/ kT. The correlation func- 
tion is complex- valued because the operators B{t) and 
Aifi) in general do not commute. While the antisym- 
metric part of CBA{t) is directly related to the response 
function by linear response theory, the power spectrum 
of the symmetrized correlation function 

SuAit) = ]^<5B{t)5A{Q) + 5A{Q)5B{t)> (4) 

depends on the Fourier transform of the dissipative part 
of the response function via 



Sba{u^) = ?icoth (^^^ X%a{^) 



(5) 



This result is the quantum version of the fluctuation- 
dissipation theorem as it relates the fluctuations de- 
scribed by Sba{^) to the dissipative part Xsa('^) of the 
response. 

In the spirit of the work by Nyquist and Johnson we 
consider as an example the response of a current 51 
through an electric circuit subject to a voltage change 
6V . This implies B — I and, because the voltage couples 
to the charge Q, A — Q. The response of the circuit 
is determined by 51{uj) = Y (ijj)5V (oj) where the admit- 
tance Y{lo) is identical to the susceptibility x/q(w). As a 
consequence of / = Q, the symmetrized power spectrum 
of the current fluctuations is given by Sjj{uj) = iujSjQ{ijj) 
so that we obtain 



Sii{u) = ?iwcoth Rer(w) 



(6) 



flLO 



huj 



ReY{Lu) 



In the high temperature limit kT ^ huj, we re- 
cover the results of Nyquist and Johnson, i.e. Sii{uj) — > 
2kTReY{u!). For the Markovian limit of an ohmic re- 
sistor, where Y{uj) = 1/R this result simplifies to read 
Sii{ll!) — 2kT/R. The quantum version was already an- 
ticipated by Nyquist in the last paragraph of his 1928 
paper—. However, he made use of the original expres- 
sion of Planck which yields only the second contribution 
present in the lower line of Nyquist thus missed the 
first term arising from the vacuum energy which already 
appears in a paper by Planck published in 1911i^ 

On the other hand, in the extreme quantum limit 
kT <C huj, we find that Sii{uj) hu}ReY{uj). In par- 
ticular, this implies that at zero frequency the spectral 
weight of the current fluctuations vanishes in the generic 
case where the admittance does not exhibit an infrared 
divergence. 

We cannot emphasize enough that the quantum 
fluctuation-dissipation relation ^ and corresponding 
implications hold true for any isolated, closed quantum 
system. Thus, upon contracting the dynamics in full 
phase space onto a reduced description of an open quan- 
tum system exhibiting dissipation these relations hold 
true nevertheless. Therefore, care must be taken when 
invoking approximations in order to avoid any violation 
of these rigorous relations. We next consider the role 
of quantum dissipation for an exactly solvable situation: 
the damped quantum harmonic oscillator dynamics. 



III. QUANTUM DISSIPATION: THE DAMPED 
HARMONIC OSCILLATOR 

A. Equilibrium correlation functions 

Let us next consider the most fundamental case of a 
simple open quantum system, namely the damped har- 
monic oscillator. This problem could be tackled by set- 
ting up a microscopic model describing the coupling to 
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environmental degrees of freedom to which energy can 
be transferred irreversibly, thus giving rise to dissipation. 
Such an approach will be introduced in Sect. IIV.AI On 
the other hand, the linearity of the damped harmonic 
oscillator allows us as alternative to proceed on a phe- 
nomcnological level. This approach is closely related to 
the usual classical procedure where damping is frequently 
introduced by adding in the equation of motion a force 
proportional to the velocity. 

Classically, the motion of a harmonic oscillator subject 
to linear friction is determined by 



Mq + M 



ds-/{t - s)q{s) + MujIq = . (7) 



In the example of an electric circuit mentioned in the 
previous section, a damping kernel j{t) with memory 
would correspond to a frequency-dependent admittance. 
In the special case of ohmic friction corresponding to 
Y{ui) = the damping force is proportional to the 

velocity of the harmonic oscillator, so that the equation 
of motion reads 



Mq + Mjq + Muolq = . 



(8) 



In Q and © the mass, frequency, and position of the 
oscillator are denoted by M, ^^nd g, respectively. Due 
to the Ehrenfest theorem, the equation of motion 10 is 
still valid in the quantum regime if we replace q by its 
expectation value. As a consequence, the quantum me- 
chanical dynamic susceptibility agrees with the classical 
expressioniSii& 



1 



1 



(9) 



where 7(0;) denotes the Fourier transform of the damping 
kernel 7(i). 

As mentioned before, the response function directly 
yields the antisymmetric part of the position autocorre- 
lation function Cqq (i) . It therefore suffices to discuss the 
symmetrized part Sqq{t) defined according to Fur- 
thermore, our linear system with linear damping repre- 
sents a stationary Gaussian process so that all higher 
order correlation functions may be expressed in terms of 
second order correlation functions^ In addition, equilib- 
rium correlation functions containing momentum opera- 
tors p can be reduced to position correlation functions 
by means oi p — Mq. The dynamics of the damped har- 
monic oscillator can therefore entirely be described by 
the response function, i.e. the Fourier transform of 

and Sqq{t)&l^ 

In the case of ohmic damping, 7(0;) — 7, the po- 
sition autocorrelation can be explicitly evaluated from 
the fluctuation-dissipation theorem |(SJ). The inverse 
Fourier transform into the time domain is determined 
by the poles on the right-hand side. The dissipative 
part of the dynamic susceptibility leads to four poles at 
uj — ±((Z> ± 17/2) with oj = (wq - 7^/4)^/^ which con- 
tribute to the correlation function Sqq(t) at all tempera- 
tures. At sufficiently low temperatures, the poles of the 



hyperbolic cotangent at w = ±w„ with the Matsubara 
frequencies Vn — 2TTn/h(3 become important as well. Af- 
ter performing the contour integration in (O , one arrives 



Sqq{t) 



2Muj 



exp 



s,mh{h(3Lu) cos(wi) -I- s\ii{h[3^ /2) sin(w|i|) 



2j_ 

Ml3^^{vl+ul) 



cosh(?i/3ijj) — coa{hPj/2) 
t'„exp(-f„|i|) 



2^,2 



(10) 



In the limit of high temperatures the second term van- 
ishes and the first term yields the classical correlation 
function. Quantum corrections to this term are relevant 
at temperatures of the order huo/k or below, and these 
corrections may be obtained from weak coupling theo- 
ries like the quantum master equation approachpiiiiSiiSiSli. 
However, there is another regime at temperatures below 
hj/ATrk. Here, the second term may initially be small, 
but nevertheless it may dominate the long-time behav- 
ior of the correlation function. This becomes particu- 
larly apparent in the limit of zero temperature where 
the exponential functions in the second term in eq. I|10|) 
sum up to an algebraic long-time behavior, i.e. Sqq{t) — 
— {fi'y /'KMuj'^)t~'^ . Its relevance for the dynamical evolu- 
tion of the damped harmonic oscillator depends on the 
details of the initial preparationSi. Although the alge- 
braic decay results from the zero temperature limit it can 
also be observed at low, but finite temperatures during 
intermediate times before an exponential decay with time 
constant vi sets \vm^. The occurrence of additional time 
scales besides 7 at low temperatures leads to shortcom- 
ings with the quantum regression hypothesis and allows 
for the decay of correlations on time scales longer than 
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B. The reduced density matrix and tlie partition function 

In the previous section, we have seen that the dynamics 
of a damped harmonic oscillator can be fully described in 
terms of the position autocorrelation function H10|l and its 
time derivatives as well as the (classical) response func- 
tion. If one is interested only in equilibrium expectation 
values of arbitrary operators acting in the Hilbert space 
of the harmonic oscillator, it is sufficient to know the re- 
duced density matrix. By means of arguments analogous 
to the dynamic case presented in the previous section, the 
reduced density matrix can only depend on second mo- 
ments of position and momentum, <q^>p and <p'^>p, 
respectively. The equilibrium density matrix then neces- 
sarily takes the form& 



1 



(27r<92>)i/2 

jq + q'? 

X cxp 



(11) 



{q 



4 



The second moments are found to read 



<q >p 



-1 -T'^ 

Mf3 ^ 



1 



and 



n M \ - 



+ I'l + Wn\l{\Vn\) 



^0 + \^n\l{\yn\) 



^1 + |i^n|7(|i^n|) ' 



(12) 



(13) 



where we have introduced the Laplace transform of the 
damping kernel 



7(^) 



(ii exp(— zi)7(t) . 



(14) 



We note that for strictly ohmic damping the second 
moment of the momentum (|13f) exhibits a logarithmic 
divergence which can be removed by introducing a fi- 
nite memory to the damping mechanism. For finite cou- 
pling to the environment, i.e. for finite damping strength 
7, the reduced density matrix (|ll|l obviously does not 
agree with the canonical density matrix exjp{—(3Hs) at 
the same temperature, where Hs denotes the Hamilto- 
nian of the undamped harmonic oscillator. 

In order to get an idea of the deviation of the true 
reduced density matrix from the canonical one, we con- 
sider the leading corrections to the second moment of 
the position due to the finite coupling to the environ- 
ment. Expanding in orders of the damping strength 7, 
we obtain for ohmic damping 



<g^>/3(7) 

<g^>/3(7 = 0) 



l + ^A, + 0(72) (15) 



with 



coth(^ 



(16) 



Here, ijj' denotes the first derivative of the digamma func- 
tion. The correction A, is depicted in Fig. ^ as a func- 
tion of temperature. We find that the leading correc- 
tions are particularly important in the quantum regime, 
kT <^ huo, while in the classical regime the corrections 
to the canonical density matrix are negligible. 

As we have already mentioned, a finite memory time 
of the damping kernel or, equivalently, a finite cutoff 
frequency cj/^ for the environmental mode spectrum is 
needed in order to keep the second moment of the mo- 
mentum (I13|l finite. If lod 3> '^Oj7i the corrections to 
the canonical density matrix for weak coupling will only 
be small if the temperature is larger than the cutoff fre- 
quency, i.e. kT ^ hujD. 

The differences between the correct reduced density 
matrix (|ll|l and the canonical one are also reffected in 
the partition function. Without specifying a microscopic 




kT/hcJo 

FIG. 1 The weak coupling correction to <q^>0 according 
to II6I I is depicted as a function of the temperature T. For 
kT 3> hujo, the correction becomes negligible. 



model for the environment, the partition function Z for 
the damped harmonic oscillator can be obtained by the 
requirement that it generates the second moment of po- 
sition according to^ 



1 



HZ) . 



(17) 



This leads to the product representation of the partition 
function, i.e. 



z^^T] 



(18) 



The properties of this partition function become more 
transparent if one relates it to a density of states p{E) 
according to2^ 



dEp{E) exp(-/3£;) 



(19) 



The factor {fi/SojQ)^^ in H18I) can then be interpreted in 
terms of the average density of states (hujo)^^ indicated 
in Fig. [3 as dotted line. We further note that the parti- 
tion function diverges for purely ohmic damping. How- 
ever, it can be shown that this divergence is entirely due 
to a divergence of the ground state energy eg in the pres- 
ence of ohmic dissipation^^. For large cutoff frequencies, 
the poles of the partition function, which determine the 
density of states, can then be determined from the condi- 



tion i/^ 



0. These poles give rise to a density 



of states which for weak damping exhibits narrow peaks 
whose width is in agreement with the result from Fermi's 
golden rule. Fig. [5] depicts an example for j/uso = 0.1. In 
view of the remark made before, the density of states is 
shifted by the ground state energy. In addition, a delta 
peak at the ground state energy has been omitted. With 
increasing damping strength, the peaks broaden so that 
for sufficiently strong damping a rather featureless den- 
sity of states results which decreases with increasin g en - 
ergy to the average density of states (cf. Fig. 3 in Ref. 23j. 
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FIG. 2 The density of states defined by inversion of the re- 
lation l|19|l for a weakly damped harmonic oscillator with 
7 = O.lojo exhibits broadened peaks close to the energies 
eo + nhjjjo. A delta function at the ground state energy eo 
is not shown explicitly. The dashed line represents the aver- 
age density of states. 



IV. DISSIPATION IN NONLINEAR QUANTUM 
SYSTEMS: THE GENERALIZED QUANTUM LANGEVIN 
EQUATION (QLE) 

A. Bath of oscillators 

For nonlinear systems the arguments given in the pre- 
vious section no longer apply. In particular, second-order 
correlation functions are not sufficient anymore to com- 
pletely describe the damped system. An alternative ap- 
proach to quantum dissipative systems starting from a 
Hamiltonian at first sight does not seem feasible because 
the absence of time-dependent forces implies energy con- 
servation. However, as we will see below, once it is real- 
ized that dissipation arises from the coupling to other de- 
grees of freedom, it is straightforward to model a damped 
quantum system in terms of a Hamiltonian. 

A well known technique to describe a statistical dy- 
namics governed by fluctuations is given by the method 
of generalized master equations and the methodology of 
generalized Langevin equations. This strategy is by now 
well developed for thermal equilibrium systems. Here 
the projector operator methodologyi2*iS*i2iSfli2& yields a 
clear-cut method to obtain the formal equations, either 
for the rate of change of the probability, i.e. the gen- 
eralized quantum master equation (QME) or the (gen- 
erally nonlinear) generalized quantum Langevin equation 
(QLE). 

Already for the case of relaxation towards a unique 
thermal equilibrium specified by a single temperature T, 
the equivalence between the two approaches is not very 
transparenlj^^. A crucial role is played by the fluctua- 
tional force which explicitly enters the equivalence, such 
as corresponding cumulant averages to an arbitrary high 
order. This fact is not appreciated generally, because one 
often restricts the discussion to the first two cumulants 
only, namely the average and its auto-correlation. It is 
a fact that little is known about the connection of the 
generalized master equation and the corresponding gen- 



eralized Langevin equation in a nonlinear situation2Si2i. 

A popular model for the dynamics of a dissipative 
quantum system subject to quantum Brownian noise is 
obtained by coupling the system of interest to a bath of 
harmonic oscillators. Accordingly, we write for the total 
Hamiltonian 



(20) 



N r 

E 



Pt 

2m, 



2 2 



qc^Xi 



where the first two terms describe the system as a particle 
of mass M moving in a generally time-dependent poten- 
tial V{q,t). The sum contains the Hamiltonian for a set 
of N harmonic oscillators which are bi-linearly coupled 
with strength q to the system. Finally, the last term, 
which depends only on the system coordinate, represents 
a potential renormalization term which is needed to en- 
sure that V{q, t) remains the bare potential. This Hamil- 
tonian has been studied since the early 60's for systems 

which are weakly coupled to the environmental degrees 
of freedom9ii8ii9i20,28i29,3Q,3i^ Q^iy j^f^-gj. ^980^ jt was re- 
alized by Caldeira and Leggett'^^ that this model is also 
applicable to strongly damped systems and may be em- 
ployed to describe, for example, dissipative tunnelling in 
solid state physics and chemical physics-^. 

One may convince oneself that the Hamiltonian H2U|) 
indeed models dissipation. Making use of the solution of 
the Heisenberg equations of motion for the external de- 
grees of freedom^ one derives a reduced system opera- 
tor equation of motion, the so-called generalized quantum 
Langevin equation^. 



Mq{t) + M f dsj{t - s)q{s) + ^^^^^ = ^{t) (21) 

J to 

with the damping kernel 



dq 



lit) = l{~t) 



N 



1 

M ^ 



nriiUj, 



■ cos(wit) 



(22) 



and the quantum Brownian force operator 
e(<) = -M-f{t-to)q{to) 

N , 

+ E "^H ^»(*o) cos(Wi[i - to]) 

H sm(wi[t - to]) 



(23) 



The generalized quantum Langevin equation (|21|) ap- 
pears first in a paper by Magalinskii^ who started from 
(|20|l in the absence of the potential renormalization term. 

The force operator (|23|l depends explicitly on the ini- 
tial conditions at time to of the bath position operators 
Xiito) and bath momenta (to)- The initial preparation 
of the total system, which fixes the statistical properties 
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of the bath operators and the system degrees of freee- 
dom, turns the force £,{t) into a random operator. Note 
that this operator depends not only on the bath proper- 
ties but as well on the initial system position qlto). To 
fully specify the reduced dynamics it is thus of impor- 
tance to specify the preparation procedure. This in turn 
then also fixes the statistical properties of the quantum 
Brownian noise. Clearly, in order to qualify as a stochas- 
tic force the random force S,{t) should not be biased; i.e. 
its average should be zero at all times. Moreover, this 
Brownian quantum noise should constitute a stationary 
process with time homogeneous correlations. 

Let us also introduce next the auxiliary random force 
ri{t), defined by 

Vit) - at) + Ml{t - to)q{to) (24) 

which only involves bath operators. In terms of this new 
random force the QLE (|21|l no longer assumes the form of 
an ordinary generalized Langevin equation: it now con- 
tains an inhomogeneous term "f(t—to)q(tQ), the initial slip 
terrr^^i^. This term is often neglected in the so-called 
"Markovian limit" when the friction kernel assumes the 
ohmic form 7(f) — > 2"fS{t). For a correlation- free prepa- 
ration, the initial total density matrix is given by the 
product pt — /3s(^o)Pbath, whcrc ps{to) is the initial sys- 
tem density matrix. The density matrix of the bath alone 
assumes canonical equilibrium, i.e. 



Pbath 



exp 




_Pl_ 

2mj 



I 2 2 



(25) 



with Af denoting a normalization constant. 

The statistical properties of the random force ri{t) then 
follow immediately: r](t) is a stationary Gaussian opera- 
tor noise obeying 



<^(^)>Pbath = 

Svvit - = 7;<v{t)v{s) + v{s)v{t)>f 



(26) 
(27) 



N 



— - — COS (uJi{t — s)) coth 



2 '■ — ' miLOi 

2—1 



2kT 



Being an operator-valued noise, its commutator docs not 
vanish 



N 



[r]{t), r]{s)] = -ihS^ — — sin (uJi{t - s)) 



4=1 



miLJi 



(28) 



Setting for the initial position operator q{to) = q^, the 
last expression in (|27|l is also valid for the noise correla- 
tion S(^^ (t) of the noise force ^ (t) provided the average is 
now taken with respect to a bath density matrix which 
contains shifted oscillators. The initial preparation of the 
bath is then given by the new density matrix pbath; 

(29) 



1 



Pbath = exp 



/3 



In some physical situations a microscopic model for the 
external degrees of freedom is available^iiM. Examples 
are the electromagnetic modes in a resonator acting as 
a reservoir or the dissipation arising from quasi-particle 
tunnelling through Josephson junctions'^^. In the case of 
an electrical circuit containing a resistor one may use the 
classical equation of motion to obtain the damping kernel 
and model the environment accordingly. This approach 
has been used e.g. to model Ohmic dissipation in Joseph- 
son junctions in order to study its influence on tunnelling 
processes^, and to describe the influence of an external 
impedance in the charge dynamics of ultrasmall tunnel 
junctions'^^. 

This scheme of the QLE can also be extended to the 
nonequilibrium case with the system attached to two 
baths of different temperature^^. Two most recent ap- 
plications address the problem of the thermal conduc- 
tance through molecular wires that are coupled to leads 
of different temperature. Then the heat current assumes 
a form similar to the Landauer formula for electronic 
transport: The heat current is given in terms of a trans- 
mission factor times the difference of corresponding Bose 
functiongSS. 

Furthermore, the QLE concept can also be extended 
to fermionic systems coupled to electron reservoirs and 
which, in addition, may be exposed to time-dependent 
driving**^. The corresponding Gaussian quantum noise is 
now composed of fermion annihilation operators. 



B. Consequences of time-reversal symmetry 

Let us now discuss some further properties of this QLE. 
If the potential V{q, t) in H2U|) does not explicitly de- 
pend on time the dynamics of the full Hamiltonian 
(|20|l obeys time reversal symmetry. It is thus an imme- 
diate consequence that the reduced dynamics must be 
invariant under time reversal as well. This must hold 
true despite the fact that the QLE has been constructed 
to allow for a description of quantum dissipation. It is 
thus instructive to see how the validity under time rever- 
sal emerges from the contracted description in terms of 
the QLE in (EB- 

Given the time of preparation t^, reversing the time 
amounts to substituting time t hy t^ — [t — to) ~ 2to — t. 
Using again the random force 77 (i) we can recast the QLE 
dynamics after the time reversal into the form 

Z"^*""* dV(q) 
Mq{2to-t)+M ds'^{2to~t-s)q{s)+ ^' 



J to dq 
- e(2<o - t) 

= 77(2to - t) - 7(2*0 - t - to)q{to) . (30) 



Ci 



■qo 



Setting next x{t) = q{2to—t) and observing that x{t) = 
-q{2tQ—t), x{t) = q(2tQ—t), we find after the substitution 
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of the integration time u = 2to — s from (|30|l the result 



Mx{t) +M f du 7(w - t)x{u) + '^^^'^^ 

J to 

= e(2io - t) 



dq 



^ rji2to - t) - j{to - t)x{to) . (31) 

Noting that the damping kernel is an even function 
of its argument, j{u — t) = j{t — u), and that a;(to) = 
q(io), we find upon changing all signs of the initial mo- 
menta Piito) —Pi{to) for the noise forces the relations 
r]{2to-t) = r]{t) and ^{2to - = ^(t). We conclude that 
the time reversed motion x{t) = q{2to — t) indeed obeys 
again a QLE of the form H21|) . This even holds true in 
the Markovian limit where j{t — s) = 2^6{t — s) as one 
can convince oneself by smearing out the delta function 
symmetrically. The QLE then reads for all times t 



Mq{t)+sgn{t-to)M-fq{t) 



dVjq) 
dq 



= m , (32) 



where sgn(x) denotes the sign of x. 

The dissipation is reflected by the fact that for times 
t > to the reduced dynamics for q{t) exhibits a damped 
(quantum)-behavior on a time scale given by the Poincare 
recurrence time^i*^; the latter reaches essentially infinity 
for all practical purposes if only the bath consists of a 
sizable number of bath oscillator degrees of freedor 
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C. Subtleties and pitfalls 

The use of the generally nonlinear QLE 121|) is limited 
in practice for several reasons. Moreover, the applica- 
tion of the QLE bears some subtleties and pitfalls which 
must be observed when making approximations. Some 
important features are: 

• The QLE (|21|l is an operator equation that acts in 
the full Hilbert space of system and bath. The cou- 
pling between system and environment also implies 
an entanglement upon time evolution even for the 
case of an initially factorizing full density matrix. 
Together with the commutator property of quan- 
tum Brownian motion, see eq. I|28|) . we find that 
the reduced, dissipative dynamics of the position 
operator q(t) and momentum operator p{t) obey 
the Heisenberg uncertainty relation for all times. 

This latter feature is crucial. For example, the non- 
Markovian (colored) Gaussian quantum noise with real- 
valued correlation S^^{t) = S^^{—t) cannot simply be 
substituted by a classical non-Markovian Gaussian noise 
force which identically obeys the correlation properties 
of (Gaussian) quantum noise S,{t). An approximation of 
this type clearly would not satisfy the commutator prop- 
erty for position and conjugate momentum of the system 
degrees of freedom. 



The literature is full of various such attempts wherein 
one approximates the quantum features by correspond- 
ing colored classical noise sources, e.g. see Refs. Il^li^. 
lia Such schemes work at best near a quasi-classical 
limi1ii^ii&, but even then care must be exercised. For 
example, for problems that exhibit an exponential sen- 
sitivity, such as the dissipative decay of a metastable 
state discussed in the next section, such an approach 
gives no exact agreement with the quantum dissipative 
theory2ii24. It is only in the classical high temperature 
limit, where the commutator structure of quantum me- 
chanics no longer influences the result. Perfect agreement 
is only achieved in the classical limit. 

The study of quantum friction in a nonlinear quan- 
tum system by means of the QLE (|21|l is plagued by 
the fact that the nonlinearity forbids an explicit solution. 
This solution, however, is needed to obtain the statistical 
properties such as mean values and correlation functions. 
This (unknown) nonlinear response function also deter- 
mines the derivation of the rate of change of the reduced 
density operator, i.e. the QME, and its solution of the 
open quantum system. 

The very fact that the QLE acts in full Hilbert space of 
system and environment also needs to be distinguished 
from the classical case of a generalized Langevin equa- 
tion. There, the stochastic dynamics acts solely on the 
state space of the system dynamics with the (classical) 
noise properties specified a prioriii. 

• The quantum noise correlations can, despite the 
explicit microscopic expression given in 127|) . be 
expressed solely by the macroscopic friction kernel 
lit)- 

This result follows upon noting that the Laplace trans- 
form 7(z) of the macroscopic friction assumes with Rez > 
the form 



7(z) 



N 



1 



2M ^ miujj 



1 



1 



z + ibJi 



(33) 



With help of the well known relation l/{x + iO+) = 
P[\/x) — iiTS{x) we find that 

Re7(z = -iiv + 0+) 

N 2 

2m '^-^ rriiUjf 

i—l ^ 

By means of (|27f) we then find the useful relation 
S^^it) = S,,{t) (35) 

/ huj 



— / dcuKe^ {— iuj + 0^) hco coth 
Jo 



\2kT 



COs{LUt) . 



In the classical limit this relation reduces, indepen- 
dent of the preparation of the bath with p or p, to the 
non-Markovian Einstein relation (t) = MkTj{t). The 
relation H35|) is by no means obvious: It implies that a 
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modelling of quantum dissipation is possible in terms of 
macroscopic quantities such as the friction kernel "/(t) 
and the temperature T. For other coupling schemes be- 
tween system and bath we generally can no longer express 
the correlation of quantum noise exclusively in terms of 
macroscopic transport coefficients. As an example we 
mention the coupling of the system to a bath of two-level 
systems (spin bath) rather than to a bath of harmonic 
oscillators^. 

Note also the following differences to the classical sit- 
uation of a generalized Langevin equation. 

• The quantum noise ^(t) is correlated with the ini- 
tial position operator q{to).'^^ This feature that 
<q{to)£,{t)>p ^ follows from the exphcit form of 
the quantum noise ^(t). The correlation function 
vanishes only in the classical limit. Note also that 
the expectation value of the system-bath interac- 
tion is finite at zero temperature. These features 
reflect the fact that at absolute zero temperature 
the coupling induces a non-vanishing decoherence 
via the zero-point fluctuations. 

Moreover, 

• the initial slip term j(t — to)'z(^o) appears also in 
the absence of the potential renormalization in the 
Hamiltonian H2U|) . With this initial value contribu- 
tion being absorbed into the quantum fluctuation 
^{t), these become stationary fluctuations with re- 
spect to the initial density operator of the bath 
Phath given by (j^ . Note, however, that with re- 
spect to an average over the bare, non-shifted bath 
density operator phath, the quantum fluctuations 
^{t) would become non-stationary. 

It is also worthwhile to point out here that this ini- 
tial value term in the QLE should not be confused 
with the initial value term that enters the correspond- 
ing QME^^'-'^*. In the case of a classical reduced dy- 
namics it is always possible - by use of a correspond- 
ing projection operator - to formally eliminate this ini- 
tial, inhomogeneous contribution in the generalized mas- 
ter equationi^'^.. This in turn renders the time evolu- 
tion of the reduced probability a truly linear dynamics. 
This property no longer holds for the reduced quantum 
dynamicE^^: For a non-factorizing initial preparation of 
system and bath this initial value contribution in the 
QME generally is finite and presents a true nonlinearity 
for the time evolution law of the open quantum dynam- 
ics! 

There exist even further subtleties which are worth- 
while to point out. The friction enters formally the 
QLE just in the same way as in the classical generalized 
Langevin equation. In particular, a time-dependent po- 
tential V{q,t) leaves this friction kernel invariant in the 
QLE. In contrast to the classical Markovian case, how- 
ever, where the friction enters the corresponding Fokker- 
Planck dynamics independent of the time scale of driving, 
this is no longer valid for the generalized quantum master 



equation dynamics of the corresponding reduced density 
matrixi^^i^^ 

The solution of the QLE involves the explicit time- 
dependence of both the friction and the potential forces. 
These in turn determine the statistical properties of the 
density matrix. As a consequence, the friction force en- 
ters the QME in a rather complex manner. This can 
already be verified explicitly for a parametric dissipative 
oscillator dynamics, where the time-dependent driving 
enters the diffusive kinetic evolution law of the reduced 
density operator or its equivalent Wigner transformSiS. 

For the bilinear system-bath interaction with the bath 
composed of harmonic oscillators it was possible to inte- 
grate out the degrees of freedom of the bath explicitly. 
Does this hold as well for other interactions? The elimi- 
nation of the bath degrees of freedom is still possible for 
a nonlinear coupling to a bath of harmonic oscillators if 
the system part of the coupling is replaced by a nonlinear 
operator-valued function of either the momentum or po- 
sition degree of freedom of the system as long as the bath 
degrees of freedom appear linearly. The resulting friction 
kernel then appears as a nonlinear friction but the influ- 
ence of the bath degrees of freedom still is obtained in 
exact formal. 

Yet another situation for which one can derive an exact 
QLE is when a nonlinear system, such as a spin degree 
of freedom, interacts with a collection of quantum (Bose) 
oscillators in such a way that the interaction Hamiltonian 
commutes with the system Hamiltonian, thus constituing 
a quantum non-demolition interaction. This case corre- 
sponds to pure dephasing and was addressed by Luczka 
for the problem of a spin in contact with a thermal heat 
bathii. It has since been rederived many times, see e.g. 
Ref.il 

We end this subsection by mentioning also the coupling 
of a system to a bath of independent fermions with in- 
finitely many excitation energies. A suitable transforma- 
tion then allows to map the dissipation onto a bosonic en- 
vironment with an appropriate coupling strengthJi^S*^. 



V. PATH INTEGRALS AND EFFECTIVE ACTION 

A. Nonlocal effective action 

A most effective approach to describe dissipation is 
based on the path integral formulation of quantum 
mechanics^^. In the path integral formulation of quan- 
tum mechanics the propagator is expressed as 

<'"l"-(,-s'")l''>^i,(0) = ,. 

(36) 

where the integral runs over all possible paths starting at 
Qi and ending after time t at qf. The paths are weighted 
with a phase factor which contains the classical action 
S[q]. 
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For the description of quantum dissipative systems it 
is important to realize the analogy between the propa- 
gator and the equilibrium density matrix. The latter is 
obtained by replacing t by —ih(3. We thus obtain from 
(|36|l the path integral representation of the equilibrium 
density matrix 



1 MTt-P) =q / 1 ^ 

5; 



where Zfs is the partition function. This integral is called 
imaginary-time path integral in contrast to the real-time 
path integral H^fifl . Note that in (|37|l the action S[q] has 
been replaced by the so-called Euclidean action S^[q] 
which is obtained by changing the sign of the poten- 
tial term as a consequence of the transition to imaginary 
times. In imaginary time we therefore have to consider 
the motion in the inverted potential. 

The connection between classical and quantum me- 
chanics becomes particularly apparent in the path inte- 
gral formulation. The dominant contribution to the inte- 
grals in (ESI and ^ arise from the stationary points of 
the action, i.e. the classical paths. Quantum effects have 
their origin in fluctuations around the classical paths. 
Therefore, it is useful to decompose a general path into 
the classical path and a fluctuation around it. Expanding 
the action in powers of the fluctuations the second order 
term yields the leading quantum corrections. Higher or- 
der terms are often neglected within a semiclassical ap- 
proximation which becomes exact for linear systems. 

In the previous section we have derived an effective 
equation of motion for the system variable by eliminating 
the external degrees of freedom. The same procedure 
may of course also be carried out within the path integral 
formalisn&^ii^^i^. The influence of the environment is 
then contained in an effective action which has to be 
added to the action of the system and which in imaginary 
time is given bjiSt^LM 

where 



^(t) = ^ |i^„|7(|i^„|)exp(zj/„T) (39) 



+ 00 



n— — oo 



and 7(z) denotes the Laplace transform of the damping 
kernel j{t). The effective action (|38|l is clearly nonlocal 
and can thus not be expressed in terms of a potential. 
If the potential renormalization term in the Hamiltonian 
(Pn|l would be absent, there would have been a local con- 
tribution in (|38() . The self interaction of the paths in- 
duced by Ip!^ via the kernel decays for ohmic damp- 
ing only algebraically as and therefore represents a 
long range interaction. 




FIG. 3 Cubic potential as defined in Eq. lHUl . 



B. Application: The dissipative decay of a metastable state 

A local potential minimum may be metastable due to 
the environmental coupling and quantum effects. Cor- 
respondingly, there are two escape mechanisms: ther- 
mal activation which dominates at high temperatures 
and quantum tunnelling which becomes important at low 
temperatures. To be definite, we consider the cubic po- 
tential 



Viq) 



M 



2 2 



1 - 



(40) 



which is depicted in Fig. 13 The barrier height is given by 
Vb = {2/27)MujgqQ and, in this special case, the barrier 
angular frequency uj^ equals the well angular frequency 
ujo. 

In Fig.^ the decay rate is shown in an Arrhenius plot. 
At the so-called crossover temperature Tq, see Eq. 
below, there is a rather distinct transition between the 
thermal regime on the left side and the quantum regime 
on the right side^. Furthermore, we observe that the 
thermal regime is larger for stronger damping, i.e. the 
system becomes more classical. 

While a real time approach to dissipative decay is 
feasibleSLSiSiSi, a simpler alternative is provided by an 
imaginary time calculation where the partition function 
Zp is considered. Since the potential (|40|l is not bounded 
from below, it is no surprise that strictly speaking Zfj 
does not exist. From the path integral point of view 
there exists an unstable fluctuation mode around the bar- 
rier which leads to a saddle point in function space. One 
can circumvent this difficulty by performing the integra- 
tion in the direction of steepest descent. The partition 
function and as a consequence also the free energy then 
acquire an imaginary part which may be related to the 
decay rate^i*^. For details of this relation we refer the 
reader to the discussion in Ref. |6^ 

The transition between thermal and quantum regime 
can be well understood within the path integral picture 
by considering the possible classical paths of duration hp 
in the inverted cubic potential. For high temperatures or 
short imaginary times h(3 the only classical solutions are 
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FIG. 4 Arrhenius plot for the decay rate of a metastable state. 
The damping strength varies from the upper to the lower 
curve as 7/2a;o = 0, 0.5, and 1 (data taken from Ref . Wi) . 



the constant solutions q ~ in the well and % = 2go/3 
at the barrier. Below a temperature given by the positive 
solution of 



t^i + |i^i|7(ki|) - = 



(41) 



a second fluctuation mode becomes unstable, thereby in- 
dicating a new classical solution which performs an os- 
cillation around the barrier—. This new solution is asso- 
ciated with quantum tunnelling. Therefore, (|41() defines 
the crossover temperature which for ohmic damping is 
given hy^ 



h 



1/2 



(42) 



As discussed above, stronger damping leads to a lower 
crossover temperature and smaller quantum regime. It 
thus makes the system more classical. A distinct feature 
of the dissipative quantum decay in the low temperature 
regime is its algebraic enhancement of the decay rate with 
temperature^^. For the case of an ohmic environment 
with a constant friction behavior at low frequencies one 
finds a universal T^-enhancement of both, the prefactor 
and the effective action, with the latter dominating the 
exponential rate enhancemenliSiiSi. 



VI. SUNDRY REMARKS AND CONCLUSIONS 

With this work we elucidated the topic of quantum 
Brownian noise which drives the dynamics of open dissi- 
pative quantum systems. We have emphasized the strong 
implications that thermal equilibrium and time-reversal 
symmetry (leading to detailed balance symmetry) im- 
poses on the reduced system dynamics. We also pointed 
out the advantageous use of the path integral scheme for 
the case of nonlinearity and strong friction. 



This method seemingly is superior to any perturba- 
tive scheme that treats the system-bath coupling to low 
orders only, such as the weak coupling master equation 
methodologjiiii2ii2i2S. There are recent developments in 
the strong friction regime, where an alternative descrip- 
tion in terms of a quantum Smoluchowski equation is 
promising^^, see also the contribution by Grabert, Anker- 
hold and Pechukas in this special issue. 

A consequent use of the so-called rotating-wave ap- 
proximations also may entail some danger. It safely can 
be applied only in the weak coupling regime for resonant 
situations. We remark that the use of the rotating-wave 
approximation implies a violation of the Ehrenfest the- 
orem in the order of 'y^^is^eg -^^jjich is clearly small only 
in the weak coupling regime, i.e. for 7 <C wq , with uJq^ 
denoting some typical time scale of the system dynamics. 
The same remarks apply to the failure of the quantum re- 
gression theoremi^iSSiii: Again, the effect might be small 
for (i) very weak damping, (ii) not too low temperatures 
obeying kT ^ and (iii) not too short evolution times. 

The generalized quantum Langevin equation discussed 
in Sec. IIVI is formally exact for nonlinear quantum sys- 
tems. Its practical use is typically restricted, however, to 
linear systems for which the response can be evaluated 
in closed form. This holds true even for time-dependent 
linear systems for which the response is still linear al- 
though the evaluation involves the use of numerical Flo- 
quet theory^. The lack of knowledge of this generally 
nonlinear response function also plagues the evaluation 
of the corresponding generalized master equations. 

This problem of obtaining the generalized master equa- 
tion from the nonlinear generalized Langevin equation 
is not solved either for the classical problem with col- 
ored noisc^^. It is also this very problem that limits the 
practical use of the various variants of recently derived 
stochastic Schrodinger approachesiS. 

Likewise, the use of nonlinear, but non-stochastic 
Schrodinger equations of the type discussed and surveyed 
in Refs. 72 73 74 can clearly not describe the time evo- 
lution of a quantum mechanical mixture, nor do these 
nonlinear deterministic approaches obey, in general, the 
Heiscnberg uncertainty relation£2i24. 

There have been repeated attempts since the early days 
of quantum mechanics to explain quantum phenomena in 
terms of Einstein's theory of classical diffusion. Early ef- 
forts in this direction were those of Fiirtb/-, Fenyes^^, 
Weizel^^ and Favella^®. This credo has been popularized 
later by Nelson^"'^^ under the label of "Stochastic Me- 
chanics" . It can convincingly be demonstrated, however, 
that a quantum dynamics is quite distinct from a clas- 
sical Markovian - or even non-Markovian - stochastic 
dynamicsiiiSS. This holds even more so, if one attempts 
to incorporate the quantum dissipation for an open sys- 
tem. 

These sundry remarks thus give clear evidence that the 
topic of quantum Brownian motion - although 100 years 
have passed since Einstein's cornerstone contributioni- 
cannot be considered as "solved". For example, little is 
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presently known also for the description and the role of 
quantum noise acting in steady state, far from equilib- 
rium situations, i.e. when several baths of different nature 
and/or different temperature are coupled to the nonlin- 
ear system of interest. 

The latter case is also of salient importance for the 
description of the quantum dynamics of so-called quan- 
tum Brownian motorsSi. In those applications quan- 
tum Brownian noise is utilized in combination with non- 
equilibrium (classical or quantum) fluctuations to per- 
form exploitable work against external bias forces. In 
summary, the field of quantum Brownian motion is very 
much alive and lots of challenges still need to be ad- 
dressed and mastered. 
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